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ON THE FLOW OF NON AXISYMMETRIC PERTURBATIONS 
OF CYLINDERS VIA SURFACE DIFFUSION 

JEREMY LECRONE AND GIERI SIMONETT 


Abstract. We study the surface diffusion flow acting on a class of general 
(non-axisymmetric) perturbations of cylinders Cr in IR^. Using tools from 
parabolic theory on uniformly regular manifolds, and maximal regularity, we 
establish existence and uniqueness of solutions to surface diffusion flow starting 
from (spatially-unbounded) surfaces defined over Cr via scalar height functions 
which are uniformly bounded away from the central cylindrical axis. Addition¬ 
ally, we show that Cr is normally stable with respect to 27r-axially-periodic 
perturbations if the radius r > l,and unstable if 0 < r < 1. Stability is also 
shown to hold in settings with axial Neumann boundary conditions. 


1. Introduction 


The surface diffusion flow is a geometric evolution law in which the normal 
velocity of a moving surface equals the Laplace-Beltrami operator of the mean 
curvature. More precisely, we assume in the following that Tq is a closed embedded 
surface in IR^. Then the surface diffusion flow is governed by the law 

(1.1) V(i) = Ar(t)Rr(t), r(0) = Tq. 

Here T = {r(i) : i > 0} is a family of hypersurfaces, V{t) denotes the veloc¬ 
ity in the normal direction of T at time t, while Ar(t) and Rr(t) stand for the 
Laplace-Beltrami operator and the mean curvature of T(t) (the sum of the princi¬ 
pal curvatures in our case), respectively. Both the normal velocity and the mean 
curvature depend on the local choice of the orientation. Here we consider the case 
where r(t) is embedded and encloses a region H(t), and we then choose the outward 
orientation, so that V is positive if H(t) grows and Rr(t) is positive if r(t) is convex 
with respect to r2(t). 


The unknown quantity in (1.1) is the position and the geometry of the surface 


r(t), which evolves in time. Hidden in the formulation of the evolution law (1.1) 
is a nonlinear partial differential equation of fourth order. This will become more 
apparent below, where the equation is stated more explicitly. 

It is an interesting and significant fact that the surface diffusion flow evolves 
surfaces in such a way that the volume enclosed by r(t) is preserved, while the 
surface area decreases (provided, of course, these quantities are finite). This follows 
from the well-known relationships 

(1.2) ^Vol(t)=/ V{t)da= [ 

Jr(t) Jt 


/r(t) 


^r(t)^r(t) da — 0 
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and 


(1.3) 


d 

dt 


A(t) — [ y{t)Hr(t)da — ( (Ar(t)i?r(t)) ^r(i) 

JT(t) JT{t) 


/r(t) 


gradr(t)77r(i) 


r(i) 


da < 0, 


where Vol(t) denotes the volume of n(t) and A{t) the surface area of T{t), respec¬ 
tively. Hence, the preferred ultimate states for the surface diffusion flow in the 
absence of geometric constraints are spheres. It is also interesting to note that the 
surface diffusion flow can be viewed as the ^-gradient flow of the area functional, 
a fact that was first observed in [5D]. This particular structure has been exploited 
in |29[ I30j for devising numerical simulations. 

The mathematical equations modeling surface diffusion go back to a paper by 
Mullins [32j from the 1950s, who was in turn motivated by earlier work of Herring 


Since then, the surface diffusion flow (1.1) has received wide attention in the 


mathematical community (and also by scientists in other fields), see for instance 
Asai [5], Baras, Duchon, and Robert [S], Bernoff, Bertozzi, and Witelski [3], Cahn 
and Taylor m. Cahn, Elliott, and Novick-Cohen [10], Davf and Curtin Elliott 
and Garcke im, Escher, Simonett, and Mayer [TS|, Escher and Mucha m, Koch 
and Lamm |25j . LeCrone and Simonett HZ], McCoy, Wheeler, and Williams m, 
and Wheeler BHl 135] . 


The primary focus of this article is the surface diffusion flow (1.1) starting from 
initial surfaces Tq which are perturbations of an infinite cylinder 

Cr ■= {(a;, y, z) € IR^ : y'^ + z'^ = r'^ , x e IR} 


with radius r > 0, centered about the cc-axis. 
The main results of this paper address 


(a) existence, uniqueness, and regularity of solutions for (1.1) for initial surfaces 
To that are (sufficiently smooth) perturbations of Cr; 


(b) existence, uniqueness, and regularity of solutions for (1.1) for initial surfaces 
To that are (sufficiently smooth) perturbations of Cr subject to periodic or 
Neumann-type boundary conditions in axial direction; 

(c) nonlinear stability and instability results for Cr with respect to perturbations 
subject to periodic or Neumann-type boundary conditions in axial direction. 

In section 2 we derive an explicit formulation of (1.1) for surfaces which are pa¬ 
rameterized over Cr in normal direction by height functions p : Cr ^ IR, see (2.5). 


Equation (2.5) is a fourth-order, quasilinear, parabolic partial differential equation 


for which well-posedness in the unbounded setting can be established by virtue of 
recent results for maximal regularity on uniformly regular Riemannian manifolds 
(c.f. [36112]). A proof of well-posedness is carried out in section 3 of the current arti¬ 
cle. Eor other recent results regarding geometric evolution equations in unbounded 
settings, refer to Giga, Seki, Umeda HU HU, wherein the mean curvature flow is 
shown to close open ends of non-compact surfaces of revolution given appropri¬ 
ate decay rates at the ends of the initial surface. The current setting differs from 
that of Giga et. al. in many ways; most notably, we consider surfaces which lack 
symmetry about the cylindrical axis; we control behavior at space infinity via 
regularity bounds (so surfaces do not become radially unbounded); and we consider 
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the dynamics of surfaces which are uniformly bounded away from the cylindrical 
axis. 

In |27j , we studied the surface diffusion flow for the special case of axisymmetric 
and 27r-periodic perturbations Fq of In this case it was shown that equilibria of 


(1.1) consist exactly of cylinders, 2fc7r-periodic unduloids and nodoids. In addition, 


we established nonlinear stability and instability results for cylinders as well as 
bifurcation results, where r serves as a bifurcation parameter. More specifically, 
it was shown that cylinders Cr are stable for r > 1, unstable for r < 1, and that 
a subcritical bifurcation of equilibria occurs at r = 1, with the bifurcating branch 
consisting exactly of the family of 27r-periodic unduloids. 

In this paper we show that the stability and instability results of [27] remain true 
for non-axisymmetric perturbations subject to periodic or Neumann-type bound¬ 
ary conditions in axial direction. In particular, we show that small 27r-periodic 
perturbations Fq of a cylinder Cr with r > I exist globally and converge exponen¬ 
tially fast to another nearby cylinder 

C{y, z, f) := {(a:, y, z) : {y- y)^ + {z- zf = f‘^, x£ M}. 

It will be shown that the radius r is uniquely determined by the volume of the solid 
enclosed by Fq and bounded by the planes a; = 0 and x = 2 'k. In order to prove 


this result, we shall show that all equilibria F of (1.1) which are 27r-periodic in 
axial direction and sufficiently close to Cr (for r > 1) are cylinders C{y,z^f), with 
(y, z, r) close to (0,0, r). This result, which is interesting by itself, is obtained by a 
center-manifold argument. In order to prove the stability result, we demonstrate 
that every cylinder Cr with radius r > 1 is normally stable with respect to 27r- 
periodic, 5c^+“-regular perturbations. As discussed in section 5, these stability 
results also encompass perturbations satisfying Neumann-type boundary conditions 
in cc-direction. 


2. Surface Diffusion Flow near Cylinders 

Throughout this section, take r > 0 fixed, and let Cr C be the (unbounded) 
cylinder of radius r which is symmetric about the x-axis. Specifically, 

Cr := {(a:,rcos(0),rsin(0)) : a; G IR,0 G T}, 

where T := [0, 27r] denotes the one-dimensional torus with 0 and 27r identified and 
T is equipped with the periodic topology generated by the metric 

di:{x,y) := min{|a; - y|,27r - jx - y|}, a;, y G T. 

Further, we equip Cr with the Riemannian metric g := dx^ -I- inherited from 

IR^ via embedding. 

Remarks 2.1. Two properties of Cr are easy to see immediately: 

(a) Cr is a uniformly regular Riemannian manifold, as defined by Amann [3]. 
Indeed, this follows from [H formula (3.3), Corollary 4.3, and Theorem 3.1], noting 
that Cr is realized as the Cartesian product of the circle rS^ and IR. On a side line 
we note that the class of uniformly regular (closed) Riemannian manifolds coincides 
with the class of closed, complete manifolds of bounded geometry, see m- 

(b) Cr is an equilibrium of <0 for every radius r > 0. 
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Indeed, noting that the mean curvature Hcr = 1/?' is constant throughout the 
manifold, it follows that 

= 0 , 

so that the normal velocity = 0 and the surface remains fixed in space. 

Now, consider a scalar-valued function p : —>■ H such that 

p{p) > —r, for all p G Cr, 

and define a new surface F = r(p) over the reference manifold Cr as 

r := {p + p{p)vip) :pGCr}, 


where v denotes the unit outer normal field over Cr- We also note that the points 
p G Cr are uniquely determined via the surface coordinates {x, d) G IR x T, and so 
we also view the height function p as a function of the variables x and 6. Thus, 
abusing notation, we will interchangeably refer to p{p) and p{x,9), where, in fact, 

p{x, 9) := p((a;, r cos(0), r sin(0))), (a;, 0) G IR x T. 

Extending this slight abuse of notation, one easily identifies an explicit parametriza- 
tion for r(p) via 

(2.1) p{x, 9) := (a;, r cos(0), r sin(0)) -|- p{x, 9)(Q, cos(9), sin(d)), (a;, d) G IR x T. 


With the above construction, we have a correspondence between height functions 
p > —r and a particular class of embedded manifolds r(p) C IR^ which we refer to as 
axially-definable — namely, those manifolds which admit an axial parametrization 
(p > 0 of the form (2.1). Within this class of manifolds, the geometry, time- 
dependent evolution, and regularity are determined by the height function p itself. 
We turn to the task of explicitly expressing these relationships, and deriving the 
governing equation for surface diffusion, in the following subsection. For the purpose 
of expressing the necessary geometric quantities, we assume for the moment that 
p = p{x, 9) is smooth enough for all following calculations to work. The precise 
desired regularity of p will be addressed in detail when we discuss existence and 
uniqueness of solutions. 


2.1. Geometry of r(p). Utilizing explicit parametrizations of the form (2.1), we 
find the coefficients of the first fundamental form of r(p) as 

5ii = (1 + pI), 9i2 = 521 = PxPe, 522 = ((?' + pY + Pe)^ 

and the metric is thus given by p = pn dx'^ + 2gi2 dxd9 + 522 d9'^. Then, 

1 + 5x 


+ P? + pI 

9 - g - ^ 


9 ^^ = 9 ^^ = 


-PxPe 


5 ^^ = 


where the matrix is the inverse of [py ] and 

( 2 . 2 ) g = g{p) ■- detipij) = {r + p)^(l -f pi) + pi- 
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Likewise, the second fundamental form is expressed via the coefficients 


I 

I 

I 


_ (r + p)p^^ 

VG ’ 

TT (^ + p)pxe - PxPe 

'-- 7 ! -’ 

_ (r + p){pe 0 -{r + p)) - 2pl 


22 — 


VQ 


The mean curvature 'H(p) ^r(p) and Laplace-Beltrami operator ^T{p) are 

expressed using well-known formulas, c. f. [inilSl]- From [311 Equation (5)], 


(2.3) nip) 


Pg - (r + p)[{{r + pY + pDpxx + (1 + pDpee - ‘^PxPePxe] 

gm 


1 


and, from |34l Section 2.7], 


Ar(p) - 


VG 


(r + pf + pI 

VQ 


d^- ^de 


de 


^ + pj 


The normal velocity of the surface r(/5) = r(p(t)) is likewise 



(r + p)pt 

'A(p) = 

Throughout, we often streamline notation by omitting explicit reference to de¬ 
pendence upon the spatial or temporal variables; e.g. p = p(t) = p{x, 9) = p{t,p) = 
pit, X, 9) and r(p) = T{p{t)), etc. 



2.2. The Surface Diffusion Flow. With the formulas from the last subsection, 
we can now express the governing equation for surface diffusion of r(p) as an evolu¬ 
tion equation for the height function p alone. Thus, defining the (formal) operator 


G{p) := 


(2.4) 


{r + p) 


P? 


-h 


VW) 

1 


d,H{p) - deH{p) 




VW) 


P- deUip) - 5,H(p) 


^/W) 


we arrive at the expression 


(2.5) 


Pt{t,p) = 

P(0) = Po, 


[Gipmip), 


iOT t > 0, p G Cr, 
on Cr, 


for the surface diffusion flow for axially-definable surfaces r(p). 


Remark 2.2. The notation [G{p{t))]{p) reflects the functional-analytic framework 
we use to address equation ( |2.5| ). For each t > 0, we consider p{t) as an element 
of an appropriate Banach space of regular functions defined on C^. Then, G maps 
p(t) to another function defined over Cr which we then evaluate pointwise using the 
notation [G'(p(<))](p). 
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3. Existence and Uniqueness of Solutions 


In this section we establish existence and uniqueness of solutions for (2.5). One 
essential tool that we use throughout is the property of maximal regularity, also 
called optimal regularity. Maximal regularity has received considerable attention in 
connection with nonlinear parabolic partial differential equations, c.f. [HIIIIIIIEHI 


3.1. Maximal Regularity. Although maximal regularity can be developed in a 
more general setting, we will focus on the setting of continuous maximal regularity 
and direct the interested reader to the references [TJ H5] for a general development 
of the theory. 

Let /re (0,1], J := [0,T], for some T > 0, and let if be a (real or complex) Ba¬ 
nach space. Following the notation of [12] , we define spaces of continuous functions 
on J := J \ {0} with prescribed singularity at 0 as 

E) := e C{j, E)-.[t^ t^-^^u{t)] € BUC{j, E) and 

(3-1) lim = ol, /re ( 0 , 1 ) 

t->-o+ J 

||w||bi_,, := supti-'"||rr(i)||£;, 
teJ 

where BUC denotes the space consisting of bounded, uniformly continuous func¬ 
tions. We also define the subspace 

BUCl_^iJ,E) := {rr e C'i(j,U) : it, u e RC/C'i_;,( J, A)} , /r e (0,1) 
and we set 

BUCoiJ, E) := BUC{J, E), BUC^{J, E) := BUC^{J, E). 

If J = [0, a) for a > 0, then we set 

Ci_^(J,E) :={rre C(j,U) :rreBt/Ci_^([ 0 ,T],£;), r<supJ}, 

Cl_^{J,E) := {u e C\j,E) :u,uG Ci_^(J,if)}, /r e (0,1], 

which we equip with the natural Frechet topologies induced by BUCi-^{[0,T], E) 
and BUCl_^{[0,T], E), respectively. 

If El and Eq are a pair of Banach spaces such that Ei is continuously embedded 
in Eq, denoted Ei ^ Eq, we set 

Eo(J):=BUCi.^(J,Eo), /re (0,1], 

Ei(J) := BUCl_^{J,Eo) n BUCi_^{J, Ei), 

where Ei(J) is a Banach space with the norm 

I|u||ei(j) := supti-'"(||u(t)||i 5 „ -h ||rr(t)||i 5 i). 
teJ 

It follows that the trace operator 7 : Ei(J) —>• Eq, defined by yu := v(0), is well- 
defined and we denote by 7 E 1 the image of 7 in if 0 , which is itself a Banach space 
when equipped with the norm 

||a::||T,Ei := inf |||t||ei(j) : v e Ei(J)and 7 r; = xj. 
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Given B G £(i?i,i?o), closed as an operator on Eq, we say (Eo(J),Ei(J)) is a 
pair of maximal regularity for B and write B G A4TZf^{Ei, Eq), if 

^ ^isom (Ei(J),Eo(J) X 7 E 1 ), ^g(0, 1), 

where £isom denotes the set of bounded linear isomorphisms. In particular, B G 
MTZi_i{Ei,Eo) if and only if for every (/, wq) € Eo(J) x 7 E 1 , there exists a unique 
solution u G Ei(J) to the inhomogeneous Cauchy problem 

u{t) + Bu{t) = f{t), tGJ, 
u(0) = Mo- 


1 + 


Moreover, in the current setting, it follows that 7 E 1 = (iloj £i)° 001 i-®- the trace 
space 7 E 1 is topologically equivalent to the noted continuous interpolation spaces 
of Da Prato and Grisvard, c.f. [HliaiTslEH]. 


We shall establish well-posedness of (2.5) in the setting of little-Holder regular 
height functions p. These results seem to constitute the first existence and unique¬ 
ness results for the surface diffusion flow acting on unbounded (closed) surfaces. 
Note that we enforce minimal conditions on the ends of the initial surface, namely 
we look at surfaces which are uniformly bounded away from the axis of definition 
(i.e. the x-axis in our current setting) and satisfy minimal regularity assumptions. 

For the convenience of the reader, we include a brief definition of the little- 
Holder spaces on Cr- For a G (0,1) and U C K" open, we define bd^iU) to be the 
closure of the bounded smooth functions BC°°{U) in the topology of the bounded 
Holder functions BC‘^{U). Further, for /c G N, we define &c^+“(C/) to be the space 
of fc-times continuously differentiable functions such that the fc**'-order derivatives 
are in hc°‘{U). We then define the space of {k -I- a) little-Holder regular functions 
onC, &c'=+“( 0 , via an atlas of local charts and a subordinate localization system. 
For more details regarding function spaces on uniformly regular (and singular) 
Riemannian manifolds see [2 |3S] . 


3.2. Quasilinear Structure and Maximal Regularity. Expanding terms of 
(2.4), it is straight-forward to see that G{p) is a fourth-order quasilinear operator 
of the form 


G{p) = -A{p)p + F{p) 


^ bp{p,d^p,d‘^p) d^p 
J/3|=3.4 


E{p,d^p,d'^p), 


where /3 = {Pi,/32) G is a multi-index, |/3| := /3i - 1-/32 its length, := 

, and d^p, for fc G N, denotes the collection of all derivatives p for |/3| = k. 

This quasilinear structure plays an important role in our well-posedness results 
below, for which we must look more closely at the fine properties of the variable 
coefficient linear operator A{p). 
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Expanding terms, we have the following highest-order variable coefficients for 
the linear operator A{p)'. 


^4,0)(P) 

^(3.1)(P) 

^( 2 , 2 ) (P) 

^(1.3)(P) 

\0A){P) 


{i.r + p?+pl? 

g2 

_ 4:PxPe{{r + p)'^ + Pg) 

g2 

2{{r + pf+pl){l + pl)+^plpl 

_ '^PxPejl + pD 
g2 

a+pir 

G^ ■ 


Therefore, the principal symbol cr[^(p)] of the linear operator A satisfies 


j + fc—4 


= E lbuMp)]ip)^i^2 

j-\-k —4 

(3-3) = ^ ^((r + p)^ + pg)^l + (1 + Px)^2 ~ 2pxPeii^'^ 

> ^ ((^ + pfil + el) • 


With these quantities explicitly expressed, we will show that A satisfies the 
property of continuous maximal regularity on Cr- Well-posedness of (2.51 then 
follows by exploiting the quasilinear structure of the parabolic equation. 


Proposition 3.1 (maximal regularity). Fixe > 0, a S (0,1), and take p S [1/2,1] 
so that 4p + a ^ Z. Further, let V)i := hA^^°'{Cr) H [p > £ — r] he the set of 
admissible height functions. Then it follows that 

{A,F) e (v^,M,{bc‘^+^{Cr),bc‘^{Cr)) x 6c“(C)), 

where C‘^ denotes the space of real analytic mappings between (open subsets of) 
Banach spaces. 


Proof. The regularity of .A(-) and F{-) is a consequence of the analyticity of the 
following mappings: (which is easy to show using standard methods in nonlinear 
analysis and theory of function spaces) 

[p 1/p] : bo bc°‘{Cr), 

[p hA d>^p] : 6c'^+l'^l(0 ^ bc'^iCr), 

[(p, h) l-A ph] : bc°‘{Cr) x bc°‘{Cr) —t hc°‘{Cr), 

where Vq := bc°‘(Cr) n [p > £ — r]. Additionally, one notes that the regularity of 
A : t C{bc'^^°‘{Cr),hc°‘{Cr)) is determined by the regularity of the coefhcients 

bp : Vfj, ^ bc°‘{Cr). Then, one is left only to show that A{p) is in the denoted 
maximal regularity class for p G Note that (6c"‘+“(Cr), 6c“(Cr)) is an 
open subset of £(6c"‘+“(Cj.), bc^{Cx)) (by [TH Lemma 2.5(a)], for instance); thus the 
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regularity of A mapping into A4iy(bc^~^°‘(Cr), bc°‘(Cr) is determined by the regularity 
of A mapping into £(5c^+“(Cr), &c“(Cr)). 

To show A(p) G M^(bc‘^~^°'(Cr), bc°‘(Cr)), by [3S1 Theorem 3.6], it suffices to show 
that A{p) is uniformly strongly elliptic on Cr', i.e. we need to show that there exist 
constants 0 < ci < C 2 such that 


Cl < Re((T[A{p)]{p, ^)) < C 2 , for {p,C) G Cr x IR^, with |^| = 1 . 

However, these bounds are obvious from the expression (3.31 and the assumption 
that admissible functions p are uniformly bounded from below and above on Cr. □ 


3.3. Well—Posedness of (2.5). We now prove existence and uniqueness of solu¬ 
tions to (2.5). Due to the non-compact setting, we must take care of the behavior 
of the height function as the axial variable approaches ±oo. For this purpose, we 
concentrate on surfaces which remain uniformly bounded away from the x-axis (i.e. 
height function p > —r -|- e for p : Cr —>■ M and e > 0). Also note that functions in 
the spaces bc'^{Cr) are bounded from above by definition. 


Proposition 3.2 (existence and uniqueness). Fix e > 0, a G (0,1), and take 
p G [1/2,1] so that Ap + a ^ Z. For each initial value 

Po G Vf, := 6c'‘'"+“(Cr) n[p> e-r], 


there exists a unique maximal solution to (2.5) 

P(-,P0) G C!_^(j(po),6c“(Cr)) nCi_^(j(po),&c4+“(Cr)), 


where J{po) '■= [0,t+(po)) C 1R_|_ denotes the maximal interval of existence for 
initial data pQ. Further, it follows that 


U J{Po)^{Po} 


is open in 1R+ X V/ and the map [(t,po) c(^)Po)] defines an analytic semiflow on 

V/. Moreover, if the solution p{-, po) satisfies: 


(i) p{-,po) G UC{J{po),bc'^^+°‘{Cr)), and 

(ii) there exists S > 0 so that distbr4^+Q(c^)(p(t, po); 5b/) > ^ for all t G J(po)- 
Then it must hold that t’^(po) = oo and so p(-,po) is a global solution of ( |2.5| ). 

Remark 3.3. Condition (ii) of the proposition can also be interpreted by explicitly 
identifying the boundary of the admissible set V)i. In particular, p(-,po) remains 
bounded away from if and only if there exists some 0 < M < oo such that, for 
all t G J(po) 


(ii.a) p{t, po){p) > 1/M -I- £ — r for all p G Cr, and 
(ii.b) ||p(t,Po)||bc'‘M+o(c,) < M. 


Proof. Existence of maximal solutions is proved in the same way as Proposition 2.2 
of [27], whereas the claim for global solutions differs slightly from Proposition 2.3 of 
m due to the fact that we do not have compact embedding of little-Hdlder spaces 
over the non-compact manifold Cr. In particular, well-posedness follows from [H 
Theorems 3.1 and 4.1(c)], jSS] Proposition 2.2(c)], and Proposition 3.1 while the 
semiflow properties follow from [T2] when p < 1 and from jS] in case p = 1. □ 
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3.4. Well—Posedness: Axially Periodic Surfaces. In order to address the sta¬ 
bility of cylinders under the flow of (2.5) — for which we will use explicit information 
regarding the spectrum of the linearization of G — we restrict our setting to one 
within which spectral calculations are tractable. 

For a > 0 fixed, define the axial shift operator 

Ta '■= \x ^ X + ci\ : IR —?► IR, 

which naturally acts on functions p € bd^lCr) as 

[Tap]{x,9) = p{x + a,e). 


We define 

bCa,periCr) ■= {p & : TaP = P}, 

which we refer to as axially a-periodic little-Holder functions on Cr- It is a straight¬ 
forward exercise to see that bc'^^periCr) forms a closed subspace of bc'^{Cr), and is 
hence a Banach space. We thus consider the properties of the surface diffusion 
evolution operator G as it acts on the axially periodic spaces. 


Proposition 3.4 (G preserves periodicity). If p & n [p > —r], then 

G{p) G bc'f p^^{Cr). 

Proof. It suffices to show that Ta commutes with the nonlinear operator G, since 
p a-periodic would then imply TaG{p) = G(Tap) = G{p), as desired. The fact that 
Ta indeed commutes with G follows directly from the commutativity of Ta with the 
operations: 

[P ^ 1/p], [P ^ d^p], and [{p, h) ^ ph], 
where, as before, □ 


Remark 3.5. The fact that G preserves periodicity can also be seen directly from 
the geometric setting. In particular, recall that G was constructed by modeling the 


evolution of the surface r(p) via (l.I), which depends only upon the geometry of 


the surface r(p) itself. Thus, if p is periodic, then r(p), and all relevant geometric 
quantities on r(p) must also be periodic, and hence the evolution cannot break this 
periodicity. 


The well-posedness results for (2.5) simplify slightly when restricted to periodic 


surfaces, owing to the fact that the evolution is determined by the restriction to 
one interval of periodicity. In essence, periodic surfaces are expressed entirely in a 
compact setting. 


Proposition 3.6 (periodic well-posedness). Fix a G (0,1), a > 0 and take p G 
[1/2,1] so that 4p + a ^ Z. For each initial value 

Po G V^.per := 6c^0,t?(Cr) F [p > -r], 
there exists a unique maximal a-periodic solution to ( |2.5[ ) 

Pi;Po) G G/_^(J(po),<pe.(C.)) nGi_p(J(po),&c4+:,(C)). 

It follows that 

'D ■= IJ J(po) X {po} 

Po ^ ^fi,per 

is open in IR+ x V^^per and the map [(t,po) t p(t,Po)] defines an analytic semiflow 
on Vfi^per- Moreover, if there exists 0 < M < oo so that, for all t G Jipo), 
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a) p{t, po){p) > 1/M — r for all p G Cr, and 

b) ||p(i,po)|lb4M+o(c^) < M, 

then it must hold that ^^{pq) = oo and p{-,po) is a global solution. 


Proof. By the assumption po > and periodicity, it is clear that po is uniformly 
bounded away from the cc-axis. According to Theorem |3.2[ the surface diffusion 


flow (2.5) admits a unique solution p(-) = p(-, po) on a maximal interval of existence 
J(po). It, thus, only remains to show that p{t) is a-periodic for each t G J{po). Let 
Pa(t) := Tap(t) for t G J{po)- Then 


dtPa{t) = Tadtp{t) = TaG{p{f)) = G{Tap{t)), t G J {po), 
Pa(0) = Tap(O) = TaPo = PO, 


showing that p(-) and Pa(-) both are solutions of ( |2.5[ ) with the same initial value. 
By uniqueness, Tap(t) = p{t), implying that pit) is a-periodic for all t G J{po)- 
The global well-posedness result follows from compactness of the embedding 
bcf/p'/j.{Cr) ^ bcf pg^{Cr) and O Theorem 4.1(d)]. □ 


For p G bdf p^^{Gr) we define the area of r(p) = {p + p(p)v{p) : p G Cr} and 
the volume of the solid enclosed by r(p), respectively, in a natural way by just 
considering the part of r(p) on an interval of periodicity in the x-direction. More 
precisely, we set 


(3.4) 


Cr,a = {(x, y,z) e IR^ : -b = r^, xG [0,a]}, 

r(P, a) = {p + p{p)v(jp) : p G Cr,a}- 
Then we have the following result. 


Proposition 3.7. Fix a G (0, l),a > 0 and take p. G [1/2,1] so that Ap + a ^ Z. 
For each initial value 

Po e := H [p > -r], 

let p = p{t,pQ) be the solution of (2.5). Then the surface diffusion flow preserves 
the volume of the region enclosed by r(p(t),a) and bounded by the planes x = 0 and 
X = a. Moreover, the flow reduees the surface area of T(pit), a). 


Proof. Although the assertions follow from (1.2|-(1.3) it will be instructive to give 
an independent proof that only relies on basic computations. A short moment of 
reflection shows that volume is given by 

1 f 


(3.5) 


Vol(p(t)) = 


'[0,a]x[0,27r] 


(r + plf))^ dOdx. 


Using the short form p = p{t,x,9), and keeping in mind that p{t,x,0) is periodic 
in (x,0) G [0,a] x [0, 27r], one immediately obtains from (2.4) 

^Vol(p(t)) = J {r + p)ptdddx = 0. 

Next we note that the surface area of r(p(t), a)) is given by 


A(p(t)) = 


'[0,a]x[0.27r] 


\/G{pit)) dOdx. 
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It follows from integration by parts (using periodicity) and formulas (2.2)-(2.4) 


dt 


A{p{t)) = 



Pt dOdx 


= y (r + p)'Hpt dOdx 


^ ~ j ^ + P^s){'^xf - 2p^pe{%x'He + (1 + dOdx 

<- j ^[{r + pfiHxf + [digf] dOdx. 


Hence A{p(t)) is non-increasing and acts as a Lyapunov functional for the evolution 
in the a-periodic setting. We also conclude that a-periodic equilibria of (2.5) 
correspond exactly to (a-periodic) surfaces of constant mean curvature. □ 


Before moving on to establish dynamic properties of solutions, we state the 
following characterization for axial periodic height functions defined on Cr- This 
characterization allows explicit access to calculations involving the linearization and 
spectrum of DG{0) in the periodic setting. 

Proposition 3.8. Given p G V)j := 5c„0,+“(Cr), it follows that p satisfies the 
Fourier series representation 

p{x,e)= /5(m,n)e^™/“e“^ (x,0 )e]RxT, 

m,nGZ 

where 

p{m,n)-=^ [ 

J[0.a]xT 

Proof. This follows from the realization of p as a function over IR x T — as discussed 
in Section]^ — and classic Fourier analysis results for periodic functions in higher 
dimensions (c.f. [Mj). □ 


4. Nonlinear Stability and Instability of Cylinders Under Periodic 

Perturbations 


For the remainder of the paper we assume a = 2 tt, and consider the behavior of 
the evolution equation (2.5) acting on 27r-axial-periodic admissible height functions 
defined over Cr (i.e. po G bc^^'^eri^r) and po > —r). Specifically, we are interested 
in the nonlinear dynamics of ( |2.5[ ) in a neighborhood of p = 0 (i.e. near the 
cylinder Cr itself). Thus, we begin in the next subsection by determining the 
spectral properties of the linearized evolution operator DG{0). 


4.1. Linearization of G. Utilizing the expression ( |2.4[ ), one readily computes the 
Frechet derivative of G, which simplifies considerably upon evaluation at the con¬ 
stant function p = 0 to the operator DG{0) G C(bc 2^.per (Cr) , &cf^,per (Cr)) given 

by 

(4.1) DG{0)h = -{dl+r-^d^}{dl+r-‘^d'^+r-^}h, h G 5c2+“g^(Cr). 

This can be seen as follows, for instance. Since we already know that G is differen¬ 
tiable, the Frechet derivative DG{0) can be computed as 


DG{0)h 


d 

de 


G{eh) 

e—0 


d 

de 


^r(eh) 

£—0 


H(0)-b Ar(o)H'(0)h. 
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Here one notes that Ar(eh)'H(0) = 0 (since H(0) = = l/u is constant) while 

Ar(o) = Ac, = 52 +r-25|, 

and 'H'(O) = —(Ac,+r“^) follow from [331 formula (31)] (noting that, by orientation 
convention, mean curvature has the opposite sign in |34j 1. 

Lemma 4.1. The spectrum of DG{0) consists entirely of eigenvalues and is given 
by 

(4.2) a{DG{0)) = {— (m^ + r~^n^^ {rnf + r~'^n^ — : to, n G Z} . 

Moreover, 0 is a semi-simple eigenvalue of multiplicity 3 and 
N{DG{0)) = span{l, cos(0), sin(0)}. 

Proof. Due to compactness of the embedding hc 2 ^^g^(Cr) into &C 2 ^ pg,(Cr) we know 
that the spectrum of DG(0) consists only of eigenvalues. Taking advantage of 
Proposition |3.8| we conclude that 

o-(DG(0)) = {— (to^ + r~‘^n^) (to^ + r~‘^n^ — r“^) : to, n G Z} . 

Recalling the expression ( |4.1| ) for the linearized evolution operator, ZIG'(O) is real¬ 
ized as a Fourier multiplier acting on /i G bc\'!f^^j.{Cr) as 

DG{0)h= ^(m,n)[(TO2-fr-2n2)(TO2-tr-2n2-r-2)]e™^e“^ 

m,nGZ 

Locating the kernel of the linearization, if DG{0)h = 0, then the fact that r > 1 
necessarily implies that h(m, n) = 0 whenever to ^ 0 , and for |n| > 1 when to = 0 . 
Thus, the kernel of DG{tS) coincides with the 3-dimensional subspace spanned by 
the first order cylindrical harmonics {1, cos(0), sin(0)}. 

Similarly, if DG{0)h G N{DG{0)), it follows that h{m,n) = 0 for all to 0 
and for |n| > 1 when to = 0. Hence, N{DG{0)^) = N{DG{0)) and 0 is thus a 
semi-simple eigenvalue of DG(0) with both geometric and algebraic multiplicity 3 
and eigenspace N. □ 

Note that DG{0) C (— cxd, 0] when r > 1, which we will exploit in order to obtain 
stability. However, accounting for the presence of 0 in the spectrum DG{0) — and 
acknowledging the fact that the family of cylinders Cr are obviously not isolated 
equilibria — we turn to a center-manifold analysis for the stability analysis. 

4.2. Local Family of Cylinders. Before we can apply the main result of [53] . 
we must address the local nature of the family of equilibria surrounding the cylin¬ 
der. We show in this section that the family of equilibria surrounding p = 0 is a 
smooth 3-dimensional manifold in &C 2 )Ji“g,(Cr). Specifically, we shall show that the 
local family of equilibria coincides precisely with the manifold of two-dimensional 
cylinders near Cr, for which we provide an explicit parametrization. 

Denote by Mcyi C bc'^~^°‘{Cr) the family of height functions p = p{y,z,f) such 
that 

r(p) = G{y, z, f) := {(t, y, z) G IR^ : (p - yf {z - zf =f^, xG IR} 

i.e. r(p) is an infinite cylinder with axis of symmetry (x,y,z) and f > 0, (in 
contrast to Cr which is centered about the axis (x, 0,0) with radius r). It follows 
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that p = p{x,9) must be independent of a;, and, as depicted in Figure 1, it is clear 
that p G M.cyi must satisfy the equation 

((r + p) cos(0) - yY + ((r + p) sm{d) - z)'^ = 



Figure 1 . Cross section of C{y, z, r) defined over Cr via the height 
function p = p{x, 9). 


Solving this quadratic equation for (r + p) readily yields the formula 

p(x, 9) = cos{9)y + sm{9)z + y/f'^ — (sin(d)y — cos(0)z)^ — r. 

The mapping [(y, z, f) i-A p] is a local diffeomorphism from some neighborhood 
O C IR^ of (0,0, r) into , so that Mcyi is a 3-dimensional Banach 

manifold of equilibria containing p = 0. 


Next, we wish show that Aicyi contains all equilibria of (2.5) in a neighborhood 
of p = 0. We accomplish this by showing that M.cyi coincides locally with a stable 
center manifold for (2.5). 


Theorem 4.2. Given r > 1, (2.5) admits a unique, 3-dimensional stable eenter 


manifold in a neighborhood of p = 0. Moreover, this eenter manifold must coincide 
with Aicyi, whieh thus contains all equilibria near p = 0. 


Proof. For N = spanjl, cos(0), sin(0)} let 7r° : bc'^{Cr) —>■ iV be the orthogonal 
projection, given by 

:= {h\l)l + {h\ cos{9)) cos{9) + {h\ sin(d)) sin(0), 

where the inner product (j-) is defined for 27r-periodic, L 2 integrable functions on 
Cr- In particular, 

{h\g) := ( h{x,9)g{x,9)dxd9 

compares functions only over one interval of periodicity along the x-axis. The 
space of periodic L 2 functions over clearly contains as a subspace, 

for any cr > 0. Defining tt® := id — and S := we get the 

direct topological decomposition = N (B S which reduces the operator 

DG{0), as follows from the fact that commutes with DG{0). We conclude that 
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tt'^ coincides with the spectral projection of DG{Q) corresponding to the spectral 
set {0}. The existence of a (local) center manifold now follows from [33 Theorem 
4.1]. In particular, there exists a neighborhood U oiO & N and a mapping 

i, e C\U, for any A: e N \ {0}, 

such that ' 0 ( 0 ) = 0 , dip{0) — 0 and the graph 

:= graph(i/-) cN(BS = 6c4+“^^(C,) 

is a locally invariant manifold for the semiflow generated by ( |2.5| ), containing all 
global small solutions (and, in particular, all equilibria close to p = 0). Moreover, 
it follows that the tangent space coincides with the eigenspace N, so 

is a 3-dimensional manifold containing p = 0. 

As contains all equilibria that are close to 0, we know that A4cyi C A4°. 
Since Aicyi and A4'^ both are (smooth) manifolds of the same dimension, a stan¬ 
dard inverse function theorem argument yields Mcyi = A4° in a sufficiently small 
neighborhood of 0. But this implies that every equilibrium close to p = 0 lies on 
Aicyi, and thus is a cylinder. □ 

4.3. Stability. We are now prepared to prove the main result of this section re¬ 
garding stability and instability of the family of cylinders Cr under sufficiently small 
admissible perturbations exhibiting 27r-axial-periodicity. 


Theorem 4.3 (stability and instability of cylinders). Let a G (0,1) and take p G 
[1/2,1] so that 4p -I- a ^ Z. 

a) (Exponential stability) Letr > 1 anduj G (0,1 —r“^) be fixed. Then there ex¬ 
ist positive constants S = d(uj) and M = M{uj) such that, given an admissible 
periodic perturbation 

Po G := bclf,'^pl^(Cr) H [p > -r], 

with ||po||6c'‘>'+° < <5) the solution p(-,po) exists globally in time, and there 
exists a p G M.cyi such that 

\\p{t,po) - p||bc4^+o < Me"‘^*||po|!6c4M+c. 


b) 


Moreover, the radius r of the limiting cylinder p is uniquely determined by 
preservation of periodic volume; in particular, f is determined by the equality 
Vol{po) = Vol(p). 

(instability) Let 0 < r < 1, then p = 0 is unstable in the topology of 


be 


4fT-\-OC 

2'K^per 


{Cr 


Proof, (a) Lemma 4.1 and Theorem 4.2 show that Cr is normally stable, i.e. we 
have with A = DG{0) 

(i) near p = 0 the set £ of equilibria constitutes a C^-manifold in &C 2 )/_“g,,(Cr.) 
of dimension 3, 

(ii) the tangent space for £ at 0 is given by N{A), 

(iii) 0 is a semi-simple eigenvalue of A, i.e. N{A) 0 R{A) = &C 2 ^ pg,,(Cr), 

(iv) a(A) \ {0} C C_ = {z G C : Rez < — w} for some w > 0 

and exponential stability follows from |331 Theorem 3.1]. To determine the radius 
f of the limit p, note that Theorem 3.7 implies the equation Vol(p) = Vol(po) must 
be satisfied. However, there is a one-to-one correspondence between the enclosed 
volume Vol(p) (as defined by (3.5|) and the radius r of p, hence f is uniquely 
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determined by volume preservation. 

(b) Nonlinear instability of cylinders Cr with radius 0 < r < 1 has already been 
shown in m Theorem 5.1]. To be precise, the authors proved this result in |27] 
within the class of axisymmetric perturbations of the cylinder. However, these 
perturbations also reside in the class of axially-definable surfaces, so nonlinear 
instability of the cylinder has already been established in the current setting. □ 


Remarks 4.4. (a) We remark that the proof of the exponential stability result can 
also be based on a center manifold analysis, as in [181 Theorem 1.2]. 


(b) There is a striking difference between the situation in [18j and the current sit¬ 
uation concerning equilibria for the surface diffusion flow. While it can readily be 
inferred that (all embedded, compact, closed) equilibria of (1.1) are spheres (single 
spheres, or multiple spheres of the same radius), the existence and classification of 
equilibria in the class of axially-definable surfaces is more difficult and complex. In 
the special case of axisymmetric and periodic surfaces, equilibria are characterized 
by the Delaunay surfaces, see m- The general case for non-axisymmetric surfaces 
seems to be wide open. 


(c) Theorem 4.3 is interesting on its own, as it characterizes all axially-definable 
periodic equilibria that are close to Cr for r > 1. 


(d) We note that the radius of the cylinder C{y, z, f) to which perturbations con¬ 
verge (when r > 1) is easily determined by the fact that surface diffusion flow 
preserves enclosed volume. However, we cannot predict the location of the central 
axis (y,z). 


5. Neumann-type boundary conditions 


It is also common to consider surfaces satisfying Neumann-type boundary con¬ 
ditions in the axial direction, c.f. laiMj. We show in this section that the well- 
posedness and stability results stated above continue to hold in subspaces of bc'^[Cr) 
satisfying such boundary conditions. 


Consider the restricted cylinder Cr^a as in (3.4) and define bcl^lfir^a) (for ct > 3) to 
be the collection of height functions p : Cr,a H satisfying the boundary conditions 

0) = dxp{a-, 9) =0 and dlp{0-\-, 9) = dlp{a-, 9) = 0, 

for all 9 gT, and such that p is bc'^ regular up to the closed boundaries of Cr,a- It 
follows directly from |5] and |3] that bc°'{Cr^a) is a Banach space of functions over 
the manifold with boundary It then follows that bc'^{Cr,a) is a closed subspace 
of 6c‘^(Cr,a), and is hence a Banach space. 


Before stating results for (2.5) in the setting of Neumann boundary conditions, 
we first note that 6c^(Cr,a) is linearly isomorphic to a closed subspace of bc 2 a^peri^r) 
via the even extension operator 


p{x,9) = ip{p){x,9) ■- 


p{2ka — x,9) if a; S [{2k — l)a, 2fca], fc G Z 
p{x — 2ka, 9) if cc G [2fca, {2k l)a], fc G Z. 


Thus, all results derived above in the axial periodic setting can also be applied to 
surfaces satisfying Neumann boundary conditions via restriction. 
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Proposition 5.1. (well-posedness for Neumann boundary conditions) Fix a € 
(0,1), a > 0 and take fj, S [3/4,1] so that 4^ + a ^ Z. For each initial value 

Po e ■= 64^+“(C,a) n [p > -r], 


there exists a unique maximal solution 


P{;P0) e Cl_^{j{po),bc%{Cr,a))nCi.^{j{po),bc^+‘^{Cr,a)). 


to (2.51. It follows that 


'^■= U J{Po)x{Po} 

P0^Vyi,N 


is open in IR+ x and the map [{t,po) i—)■ p(t,po)] defines an analytic semiflow 
on V/_ 7 v- Moreover, if there exists 0 < M < oo so that, for all t G J(po), 


a) p(t,po)(p) > 1/M — r for all p G Cr^a> o,nd 

b) < M, 

then it must hold that t~^{po) = oo and p{-,po) is a global solution. 


Proof. Given po G ^be even extension po := ipiPo) is in the admissible class 

V/_per of 2a-periodic functions on Cr, and hence Proposition 3.6 gives existence, 
uniqueness, and global well-posedness of a periodic solution 

p(lPo) G 4_4J(po),6c?,,p,,(C)) nC'i_^(J(po),6c4“,,(C.)) 


to (2.5) on Cr. It remains to show that, upon restriction to Cr,a, the solution p(t, po) 
satisfies desired Neumann boundary conditions for all t G J(po)- 
Define the reflection operator R about the x = 0 plane so that 


[Rp\{x,e) := p{-x,e). 

Then the restriction pft) := p{t)\^ satisfies the desired boundary conditions if 
and only if 

Fpx — Px^ Fpxxx — Pxxx: F{/F—aPx') — F—aPx: and R{/F—aPxxx^ — F—aPxxx ; 

where T_a is the axial shift operator discussed in section [3^ 

From a geometric perspective, it is immediately clear that reflections about x = 0 
commute with the geometric calculations used to compute the evolution operator 
G, however, confirming commutativity of R and G algebraically requires a little 
finesse. In particular, note that R commutes with the operations: 

[p 1 -^ 1/p], [p 1 -^ dgp],k G N and [(p, h) i-G ph], 

while 

However, one should observe that all odd axial derivatives of p produced in the 
operator G are found in products with an even number of similar terms. For 
instance, p^, PxPxxx, and PxPx 0 are three such products found in the expansion of 
G{p). In this way, one conhrms that the operators R and G indeed commute. 
Finally, let pFi{t) := Rpit) and PR,a(i) := RT_aP{t) for t G J{po)- Then 

dtPRft) = Rdtp(f) = RG{p(t)) = G{Rp{f)), t G J{po), 

Pr(Q) = i?p(0) = i?po = Po, 
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showing that p{-) and pr{-) are both solutions of (2.5) with the same initial value. 
By uniqueness, Rp{t) = p{t), implying that p(t) satisfies Neumann boundary con¬ 
ditions at a; = 0 for all t G J(po)- Similarly, note that PR,a{t) and T-ap(t) are both 
solutions to 


dtpit) = G{p{t)), t G J(po), p(0) = T_qPo, 
so that p{t) also satisfies Neumann boundary conditions at a; = a and thus, by 
restriction, we have a solution 


p{t) := p{t)\^ 


satisfying Neumann boundary conditions on Cr,a, as desired. 


□ 


Remarks 5.2. (a) Regarding statements of stability in the setting of Neumann 
boundary conditions, we note that stability of cylinders was derived for 27r-periodic 
perturbations of the unbounded cylinder Cr- Thus, by fixing a = tt, one can derive 
stability and instability of the bounded cylinder Cr^Tr under perturbations satisfying 
Neumann boundary conditions, with stability holding for r > 1 and instability for 
0 < r < 1 as before. 

(b) We also note that the fundamental interval of periodicity has been taken to be 
fixed at 27r throughout the article, however this is only a convenience. As discussed 
in [261 Remarks 6.2], one can take perturbations of any periodicity a > 0, and derive 
analogous stability of cylinders with radius r > a/27r and instability for cylinders 
of radius 0 < r < a/2n. 
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